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ABSTRACT 

We study the coUimation of relativistic hydrodynamic jets by the pressure of an am- 
bient medium in the limit where the jet interior has lost causal contact with its sur- 
roundings. For a jet with an ultrarelativistic equation of state and external pressure 
that decreases as a power of spherical radius, p oc r~^, the jet interior will lose causal 
contact when 77 > 2. However, the outer layers of the jet gradually coUimate toward 
the jet axis as long as ?] < 4, leading to the formation of a shocked boundary layer. 
Assuming that pressure-matching across the shock front determines the shape of the 
shock, we study the resulting structure of the jet in two ways: first by assuming that the 
pressure remains constant across the shocked boundary layer and looking for solutions 
to the shock jump equations, and then by constructing self-similar boundary-layer so- 
lutions that allow for a pressure gradient across the shocked layer. We demonstrate 
that the constant-pressure solutions can be characterized by four initial parameters 
that determine the jet shape and whether the shock closes to the axis. We show that 
self-similar solutions for the boundary layer can be constructed that exhibit a mono- 
tonic decrease in pressure across the boundary layer from the contact discontinuity 
to the shock front, and that the addition of this pressure gradient in our initial model 
generally causes the shock front to move outwards, creating a thinner boundary layer 
and decreasing the tendency of the shock to close. We discuss trends based on the 
value of the pressure power-law index 77. 

Key words: galaxies: active — galaxies: jets — hydrodynamics — relativity — shock 
waves 



1 INTRODUCTION 

The idea that AGN outflows are hig hly collimate d is supported 
by ob servations (e.g. Begelman et al.| Il984i Ijors tad et al] 
l2005h , implying that confinement must occur. The cause of 
this confinement, which may occur at distances of just a 
few tens of Schwartzschild radii from the central black hole 
fjunor et al. 1999), is however not yet well-understood. It is 
generally accepted that some coUimating agent is necessary, 
but this still allows for a variety of possibilities, both exter- 
nal (e.g. pressu re confinement via an am bient static medium 
lEichler 1982, Komissarov & Fall3ll997l inertial c onfinement 
via a slow outflow or wind surrou nding the jet Komissarov! 
11994 iBromberg & LevinsonI l2007h and internal (magnetic 
confineme nt via hoop stress due to the toroidal m agnetic field 
component lBenfordlll978l iBegelman et al.lll984h . 

Pressure confinement is of particular interest because 
accretion disk winds surrounding an AGN provide an ideal 
external medium for interaction with the jet. Moreover, an- 
other coUimating agent such as magnetic hoop stress cannot 



function alone; without an ambient medium to confine the 
globally-expansive magnetic field, collimation will not occur 
ltBegelmaijll995l: iKomissarov et al.l l2007, 20091 iKomissarovl 
I2OII ). Thus pressure confinement may be relevant both on 
its own and in conjunction with other processes. 

In this context, we note that numerical simulations that 
study confinement and acceleration of the flow by MHD pro- 
cesses treat the action of th e external medium as an ap- 
plied boundary co ndition (e.g. I Komissarov et al.ll2007ll2009l: 
lKomissarovl201lh . In this work, we study the detailed physics 
of jet collimation by the external medium as a first step to- 
wards a complete treatment of both the action of the external 
medium and magnetic effects in coUimating and accelerating 
relativistic jets. 

The location and mechanism of jet collimation is interest- 
ing because it provides insight into energy dissipation within 
the jet. In a steady jet, there are two main sites where en- 
ergy dissipation is likely to occur: at the jet spine, as a re- 
sult of kink instabilities driven by the toroidal magnetic field 
(iBegelman„1998j,,Eichler,1993.) . and at the interface between 
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the jet and its environment, as a result of shear instabilities 
dMicono et aLl l2000': 'Perucho et all 1201 d: iBarkov & Baushe\j 
l201lh or collimation shocks, as we discuss here. 

In the case of a jet with an opening half-angle of 9 and 
a bulk Lorentz factor of F, the outer edge of the jet remains 
in causal contact with the jet center only if the opening angle 
obeys F^^ > sin ^ ^ 6. In this work, we assume a relativistic 
equation of state and examine the case where causal contact 
between the jet's spine and edge has been lost. Thus we model 
the jet as having a transverse structure consisting of two com- 
ponents: an inner region in which the flow is undergoing free 
expansion as it accelerates, and an outer shocked boundary 
layer region that results from the loss of causal contact within 
the jet. The geometrical shape that the jet assumes as it propa- 
gates is a direct result of its response to the collimating forces 
exerted by the ambient medium, and calculating that shape 
for the case of pressure confinement will be a major focus of 
this paper. 

Our goal is to determine the basic jet structure and ge- 
ometry under a simple set of collimation and acceleration as- 
sumptions, providing a model for the jet's steady-state con- 
figuration. This will allow more realistic initial conditions for 
simulations of instabilities and turbulence, which will in turn 
provide a model for AGN jets to which we can compare ob- 
servational signatures. 

While we specifically reference AGN jets in this work, 
our results can easily be extended to other relativistic out- 
flows, such as a gamma-ray burst in the collapsar model, col- 
li mated by the stellar envelope it breaks through, as discussed 
in lBromberg & Levinsonll2007l (hereafter BL07). 

In this paper we use the shock conditions for a hydro- 
djTiamic, relativistic jet to derive the basic structure of the 
jet. This problem has been previously studied in the case 
of a "cold" (inertia-domina ted) jet ( Kom issarov & Fall3ll997l 
iNalewaiko & Siko"rall2009b , but we now focus on the collima- 
tion behavior of a "hot" (pressure-dominated) jet. 

We initially follow a similar approach to that of BL07, 
but deviate from its methods in our treatment of entropy dis- 
tribution within the jet. Bromberg & Levinson 2009 (BL09) 
performs a similar analjrtical inspection in the limit of small 
angles; we generalize this to all angles. 

In §2 we find solutions for the jet shape using the Kompa- 
neets approximation. In §3 we examine self-similar solutions 
for the boundary layer when a pressure gradient is allowed to 
form across the layer, and then revise our solutions from §2 
to include this pressure gradient. In §4 we conclude, summa- 
rizing the results and discussing future work. 



2 KOMPANEETS APPROXIMATION 

We consider a cylindrically symmetric, ultrarelativistic jet in- 
jected into an ambient medium that has a power-law pressure 
profile. We wrap the physics of the shocked ambient medium 
into this external pressure profile and focus on the structure 
of the jet itself In this stage of the treatment, we ignore mag- 
netic fields and assume that the external pressure due to the 
ambient medium creates the sole collimating force on the jet. 

We define R and z as dimensionless parameters in cylin- 
drical coordinates describing the radial and axial distances 
as scaled by zo, the height at which the jet initially encoun- 
ters the external medium. The jet is injected with an initial 
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Figure 1 . Diagram indicating the angles and orientation of the axes 
relative to the inner shock wall. 9j is the angle the jet streamline 
makes with the 2-axis. 'i'j is the angle it makes with the shock normal, 
and Qj is the angle between the shock normal and the vertical. 



opening half-angle of 6o and impacts the wall of the ambient 
medium at the point thus denoted as {Ro — tan^o, 20 ~ 1)- 
We assume the jet is injected from a point source with steady 
flow, and streamlines are conical and characterized by the an- 
glee,. 

We further suppose that the interior of the jet is undergo- 
ing free expansion. As relativistic adiabatic expansion obeys 
pV"'''^ ~ const, it therefore exhibits a corresponding pressure 
profile of pj oc r~*, where r is the spherical radius defined 



as r oc [R 



,2^1/2 



in cylindrical coordinates. The gradual 



acceleration via adiabatic expansion of the jet interior is gov- 
erned by the relativistic Bernoulli equation, which describes 
the conversion of thermal to kinetic energy as p oc F^"* along 
streamlines (see, e.g.. Landau & Lifshitz 1959) . 

Where the jet impacts the ambient medium with a super- 
sonic normal velocity, a shocked layer will form, as indicated 
schematically in Figure [T] The layer is bounded on the inside 
by a shock front, and on the outside by a contact discontinu- 
ity. There is no mass flux across the contact discontinuity, and 
the pressure must be matched on either side of it. Adopting 
a pressure profile for the ambient medium of ps oc r^'', for 
a parameter rj, this fixes the pressure external to the jet and 
immediately inside the contact discontinuity. 

We now ad opt the Kompane ets approximation 
( KompaneetsI Il960l: see e.g. BL07 and iKomissarov & Falld 
I1997 ), treating the pressure as a function only of z within 
the shocked layer. Thus the pressure profile of the external 
medium extends across the shocked layer with no pres- 
sure gradient in the axial direction, greatly simplifying the 
problem. 

In this paper, we focus particularly on the less-explored 
case of an external pressure profile with 2 < 77 < 4 (see e.g. 
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BL07 for an example of treatment of this regime). An r; < 2 
implies that if the jet begins in causal contact it will remain 
in causal contact, suggesting that a shocked boundary layer 
would not form. For 77 = 4 the ambient pressure profile is 
equivalent to that of free expansion, matching the pressure 
within the jet, and an ?? > 4 would resu lt in a rarefaction 
at the jet boundary (Begelma n et al.|[l98 4'). Thus the regime 
between these two values is a logical place to examine. 

Physically, this pressure profile range could describe mul- 
tiple scenarios for the confining medium. The ram pressure 
of a head-on wind decreases as p oc r~^, an d any ob liquity 
would serve to steepen that pressure profile iEichlei|[l982h . 
The range is similarly relevant if the confining medium were 
an accretion flow such as Bondi accretion [p oc r^^'^^), and 
it would not be an unrealistic range for a disk corona, or a 
stellar envelope in a GRB coUapsar model. 

The flow parameters across the inside boundary of the 
shocked layer are governed by the relativistic oblique shock- 
jump conditions. These establish conservation of mass, mo- 
mentum, energy, and tangential velocity: 



i^j ^j,x — /5s ^ 



J'- 
2 



Ws'^IPXx + Pa 



(1) 

(2) 
(3) 
(4) 



where the a;-direction is chosen perpendicular to the shock 
front and the y-direction is tangential (see Figure The sub- 
script j denotes quantities within the freely-expanding inner 
jet region, and the subscript s denotes quantities within the 
shocked boundary layer Here j3 — v/c, and w = e + p where 
e is the total proper energy density, given by e = p + 3p with p 
defined as the proper rest mass density. Considering the case 
of an ultrarelativistic gas in the regime where the jet is still 
accelerating, we assume that p <^p and the equation of state 
is p = e/3, such that w ^ 4p. 

Examining the geometry of the problem, we recognize 
that the angle of impact of the streamline with the shock 
front, sin>I'j, can be represented in terms of 6j and Uj, re- 
spectively the angle the streamline makes with the z-axis 
and the angle that the shock tangent makes with the 2;-axis. 
By noting that tanoj — dRj/dz, where Rj is the shape of 
the shock front, and assuming conical streamlines such that 
tan^j = Rj/z, we can combine the geometry of the problem 
with the general shock jump conditions to obtain a differen- 
tial equation governing the shape of the inner shock wall. 



+ l+(^) 



3^ + 1 

Pi 



(5) 



This can be solved analytically for dRj/dz and numerically 
integrated to find Rj{z). 

To solve for the shape of the contact discontinuity sepa- 
rating the shocked layer from the ambient medium, denoted 
as Rc, we assume energy-momentum conservation through 
a volume of the shocked layer, as described in BL07. As in 
BL07, we assume that the flow parameters within the shock 
depend only on the vertical distance z in order to make the 
problem analjrtically tractable. Making the further assumption 
that |/3j I ^ 1, we obtain a differential equation similar to that 
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Figure 2. Two examples of the solutions for the shock front and asso- 
ciated contact discontinuity. In the first plot, -q = 7/3; in the second 
Tj = 11/3. For both plots the remaining parameters are held constant 

atPs,o/Pj,o = 1, r^.o = 50, and 6*0 = 1/50. 



in BL07 governing the contact discontinuity: 



d ( 

TzV 



(6) 



We insert into this the known expressions for ps, Pj and 
Tj from the previous paragraphs. While a shock exists, Vs and 
j3s,z are obtained from the original shock jump conditions; if 
the shock closes to the axis then Fs and j3s,z are obtained 
by assuming adiabatic expansion within the fully-shocked jet. 
We then numerically integrate Eq ([6| simultaneously with the 
differential equation for the inner shock wall, solving for both 
Rj {z) and Rc{z). Thus we obtain the shapes of the shock front 
and the contact discontinuity in terms of the initial pressure 
ratio Ps,o/pj,o, the initial opening angle ^0, the initial Lorentz 
factor Fj.o, and the pressure power-law index r/. 

Examining the effects of varying these four parameters, 
we see that under certain conditions the external pressure can 
drive the shock front back to the jet axis, resulting in a fully- 
shocked jet. The jet is more likely to close when Ps,o/Pj,o 
is large and rj is small. Initial under- or over-pressurization 
of the jet strongly affects whether or not it closes, but for 
physical scenarios we would expect that the pressure is ap- 
proximately balanced where the jet first impacts the wall, 
Ps,o/Pj,a ~ 1. Increasing 77 can change whether or not the 
shock will reach the axis and, for cases where the jet does 
close, drives the point at which this occurs down the z-axis, 
further from the source. Two examples are shown in Figure 
[2j one in which the shock converges to the axis and one in 
which it doesn't. This demonstrates the effect that 77 can have 
when all other parameters are held constant. 

Whether or not a jet closes and, if it does, the value of 
z at the point where the shock meets the axis is dependent 
upon the value of the product ^oFj.o, rather than on 9o or 
Fj,o individually. Assuming that the initial pressure ratio and 
T] are held constant, we find that all jets with the same value of 
^oFj,o close at the same point on the ^-axis. Increasing ^oF^.o 
causes the jet to close further down the axis from the source, 
until the point where it no longer closes. This dependence 
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9 oTjfi 

Figure 3. Dependence of the jet closure point on 9oT^j,o - The different 
curves correspond to various values of ri, ranging from r/ = 13/6 on 
the right to = 3 on the left, in increments of 1/6. Here again, 
Psfi/Pjfi = 1- The closure point does not depend on the value of 
Vjfi individually. 



the streamlines as the jet colHmates would go hand in hand 
with a force inwards along the radius of curvature. Because of 
the special-relativistic length contraction along the direction 
of motion, the curvature appears to the fluid to be more ex- 
treme than it is in the lab frame by a factor of F, resulting in a 
sizable centripetal force for even slight curvature. Due to this 
centripetal force, a pressure gradient should then form across 
the boundary layer, with higher pressure at the outer wall of 
the jet and lower pressure at the shock front. 

With this in mind, we now refine the model in §2 by in- 
cluding effects of a transverse pressure gradient within the 
boundary layer. Even assuming a steady state and axial sym- 
metry, this problem intrinsically involves the solution of par- 
tial differential equations in two dimensions, so we simplily 
the problem by seeking self-similar solutions. In order to do 
so, we treat the streamlines as being very nearly conical (an 
assumption justified by the results of §2) and calculate their 
deviation from conical as a function of position within the 
boundary layer, thus examining the effects of a pressure gra- 
dient across the boundary layer. 



of the closure point on OoTj^ is illustrated in Figure [3] for 
2 < r; < 3, assuming ps,o/pj,o = 1. For > 3, ^oFj.o < 1 is 
required to produce a jet that closes. 

We find that the shape of the contact discontinuity is 
initially dependent upon the values of the four parameters 
Ps,o/Pj,o, 9(1, rj,o> and 77. For large z, however, the contact 
discontinuity asymptotes to one of two shapes: if the shock 
closes to the axis, then the flow is governed by free expan- 
sion beyond that point and the contact discontinuity takes the 
shape fle PC z^^'^, as found in prev ious works such as BL07, 
BL09 and^ Levinson & Eichlerl20 0^. If the shock does not close 
to the axis, however, then the contact discontinuity takes the 
shape Rc oc 2, in contrast to these previous works. 

The discrepancy between our work and previous studies 
appears to arise as a result of entropy treatment: in our work, 
we assume that as long as the shock has not yet closed to the 
axis, the continued addition of material into the boundary 
layer ensures that the boundary layer cannot have constant 
entropy throughout, and thus the layer is not governed by adi- 
abatic expansion. We instead solve for all quantities directly 
from the shock jump equations, without making assumptions 
about the behavior of material within the boundary layer. 

This difference in the behavior while the jet remains open 
is crucial, since this model for the jet can only be accurately 
applied in the regime in which it remains open. After the 
jet has closed, complex effects such as rarefaction waves or 
oblique shock reflections (e.g. Gomez et ak 1995J will likely 
arise, and this simple model no longer adequately describes 
the jet's behavior. 

The conical asymptote that we find in this approximation 
will be an important factor in our ability to refine this model, 
as shown in §3. 



3 SELF-SIMILAR TREATMENT OF THE BOUNDARY 
LAYER 

The Kompaneets assumption of constant pressure across the 
shocked layer is inconsistent with basic intuition: in a real- 
istic large-scale jet, one would expect that the curvature of 



3.1 Self-Similar Solutions 

We assume that the opening angle of the jet is much greater 
than 1 /r, such that causal contact has been lost. We further 
suppose that the boundary layer that forms has a thickness 
that is of order AS ~ 1 /F, ensuring that the boundary layer 
is very thin compared to the width of the jet. The radius of 
curvature of the jet is then much larger than the width of the 
boundary layer, allowing us to treat the curvature as a small 
effect. 

Given the assumption of nearly conical streamlines, we 
can in this case treat the entropy and the Bernoulli constant 
as being the same on all streamlines, implying potential flow. 
These assumptions are physically realizable if the majority 
of the material within the shock enters at approximately the 
same point near the base of the jet; we will later demonstrate 
that this assumption is self-consistent here. 

The flow within the boundary layer is governed by the 
energy equation, momentum conservation along streamlines, 
and mass conservation: 

p = Ap'^* (7) 

V-(p/3r) = (9) 

where A and B are constants. We henceforth treat the prob- 
lem in the spherical polar coordinates r, 6 and (f> for conve- 
nience in describing a jet with approximately radial stream- 
lines. As the maximum transverse speed that can be achieved 
without a shock forming is of order l/F, we can therefore 
assume that fig is of this order, and fir is of order one. Adopt- 
ing this characteristic scale, we state that ^ ~ fJ^. Writing 
out I3r and employing the fact that /3| + F^^ <C 1, we have 
I3r~l - + F-2). Using this, we then combine the flow 
equations, retaining terms only to lowest order. 

If we now assume that the external pressure is a power 
law in spherical radius r, pe oc r"'', then we can find self- 
similar solutions for the flow structure near the wall of the 
jet in the specific case of a pressure-dominated jet, or near 
enough to the jet source that the Lorentz factor is much less 
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than its asymptotic value (F <C Too)- In this limit, the equa- 
tions reduce to a pair of coupled partial differential equations 
for the Lorentz factor F and the transverse velocity Pe within 
the boundary layer: 

(10) 



r dr \V2 ) dd \V2 



d , „ . „ dBe 1 d 



We now seek self-similar solutions of the form 



-'7/4, 



v/ii 



(11) 



(12) 



where we have chosen the constant in the Bernoulli equation 
to be unity (i.e. pF'^ — 1) for simplicity. In these solutions 
g and h are functions of a similarity variable ^ that describes 
the distance from the contact discontinuity, normalized by the 
expected scale of the boundary layer, ^ oc (Sc — 6)/A9. The 
angular thickness of the boundary layer is expected to scale 
as AO = l/Tc, such that ^ oc r^^^iOc - 9), where Oc = Odr) 
is the location of the contact discontinuity. 

The fact that the streamlines at 6c must be parallel to the 
contact discontinuity, requiring that /3g{dc) = rdOc/dr, yields 
the further constraint that 



dr 



= /i(0)r-(^+''/^\ 



(13) 



Choosing the proportionality constant such that ^ is de- 
fined as 



1 



(14) 



absorbs the boundary condition into the similarity variable 
and ensures collimating solutions (such that ft(0) < 0). 

The boundary condition g{0) = 1 is enforced so that the 
pressure is matched at the contact discontinuity, but /i(0) is 
allowed to range. Assuming solutions of this form, we define 



and : 



(9M\ 
\h{Q)) 



(15) 



with ui(0) — g(0) = 1. Further defining a:: = (1 — we 
obtain a set of coupled linear ordinary differential equations 
for w{x) and q{x) that can be cast to reflect the existence of a 
critical point: 



— 2{x — q)^ 



{nw"^ + q{x - q)) 



(16) 



(17) 



fiw^ — 2{x — qY 

The critical point occurs where the denominator of these 
equations goes to zero, forcing the numerators to also go to 
zero at this point in order to prevent q and (tiw^)' from di- 
verging. 

The critical point can be physically understood as a type 
of sonic point, in analog y with the critical points discussed 
in Blandford & Payne 1982 (BP82). Setting the denominator 
in Eqs ( |16| ) and l |17| l equal to zero yields the constraint that 
r^/3| = 1/2 at the critical point, where /3e = /3- V^/l V5|. This 
constraint is mathematically equivalent to two conditions: 
that the speed of sound waves measured in the lab frame in 
the direction of must vanish (analogous to BP82's con- 
dition that the sound speed normal to the similarity surfaces 



is equal and opposite to the component of flow speed in the 
same direction), and that the sound speed along the similar- 
ity surfaces is equal to the fluid speed and carries no signals 
in this direction (analogous to BP82's condition that the wave 
signals are normal to the similarity surfaces). 

The physical solutions for q and w are those that pass 
through the critical point; the solutions that don't pass 
through the critical point display unphysical behavior, such 
as crashing and becoming double-valued. 

Requiring finite q' and {^.uP')' at the critical point results 
in three possible sets of relations among q,w, fi and x specifi- 
cally at the critical point. Each one of these sets of relations at 
the critical point applies uniquely to one of the three regions 
?7 < 8/3, ?7 = 8/3, and t] > 8/3 to produce physical solutions. 
Thus, while Eqs ( |16| | and ( |17D must generally be solved nu- 
merically, we can use these constraints at the critical point to 
help identify the solutions of interest. 

One such family of solutions occurs when fi = 1/2 and 
q{^) = w{^), i.e., when h{^) — —\/2g{^) everjwhere. The 
solutions forp oc g{iY this case are shown in Figure [4(a)| 

The most interesting result evident is that the structure 
of the solutions is divided based on the value of 77. For 77 = 
8/3 the solution is analjrtic and linear: g = w = 1 — 2^/3. 
This implies that the pressure decreases monotonically from 
the outer wall of the boundary layer (^ = 0) inward, with 
P <7(C)^ oc (1 — 2^/3)'', while the Lorentz factor F increases 
linearly inward. 

For r; > 8/3, within this special family, the solutions de- 
crease monotonically and asymptote to zero, with the steep- 
est decrease in q and w occurring for rj ^ 8/3. The curve for 
77 = 4 is roughly constant ax. q = w ^ 1. 

Solutions for r; < 8/3 within this special family are also 
possible, but they are not single-valued: at the critical point 
the derivatives diverge and the solutions do not continue 
to larger values of ^. For these solutions to be physical, the 
boundary layer would have to be truncated at a point before 
where the solutions crash and become double-valued. 

By looking beyond this special family of solutions, it is 
possible to identify a solution for every 2 < 77 < 4 wherein the 
pressure is a monotonically decreasing, single-valued func- 
tion of 5. 

In the case of 77 > 8/3, the unique physical solutions are 
the solutions that belong to the special family where = 1/2 
and q — w. For r/ < 8/3, however, the value of fi is not fixed. 
Instead, for each value of r/ there is a single value of fi that 
yields the physical solution that traverses the critical point. 
For these solutions, fi varies from /i = 1/2 for 77 = 8/3 to /i = 
00 for 77 = 2 (see Figure [5]l. When 77 < 8/3 the solutions for 
w asymptote to zero, but the solutions for q become negative 
at finite ^, as in the 77 = 8/3 case (where w and q reach zero 
simultaneously). The steepest decrease in q thus occurs for 
77 ^ 2, whereas the steepest decrease in w now occurs for 77 ^ 
8/3. The physical pressure solutions for r; < 8/3 are plotted 
in Figure [4(b)l and the physical solutions for q{^) (describing 
the spatial dependence of the transverse velocity) are plotted 
in Figure [6l 

Thus we find solutions that all display monotonically de- 
creasing pressure from the contact discontinuity across the 
boundary layer, as we would expect for a pressure gradient 
arising from centripetal acceleration. Solutions for 77 ^ 2 have 
a nearly constant pressure across the boundary layer, but the 
pressure profiles become steeper as rj increases to 8/3. The 
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Figure 4. Pressure within tiie shock as scaled by the external pressure. 
The contact discontinuity is located at ^ = 0. (a) Curves for 2 < 
r; < 4 in increments of 1/6 (beginning at rj = 13/6 and ending at 
r} = 23/6), for the special case that h{^) = —\/2g(£) everywhere. 
The value of rj increases when viewing the curves from lower left of 
the graph to the upper right. The curves for r; > 8/3 are physical, 
while the curves for < 8/3 are not unless truncated, (b) Physical 
solutions for 2 < r; < 8/3 only, when the relation between g and h 
is allowed to vary. The value of r] decreases when viewing the curves 
from lower left of the graph to the upper right. In both figures, the 
solid curve corresponds to = 8/3; the orange region contains curves 
for r] < 8/3 and the blue region contains curves for r] > 8/3. 
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Figure 5. Some values of /i as a function of rj for 2 < < 8/3, found 
empirically by seeking, for each value of rj, the solution that passes 
through the critical point. A numerical fit to these data is shown on 
the left, and a table of values is given to the right. 
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Figure 6. Physical solutions for q{^) within the shock, where the re- 
lation between g and h has been allowed to vary. The contact discon- 
tinuity is located at ^ = 0. Curves for 2 < < 4, in increments of 1 /6 
(again beginning atr} = 13/6 and ending at = 23/6), are plotted. 
The value of rj increases when viewing the curves from lower left of 
the graph to the upper right, with the solid curve corresponding to 
») = 8/3. The coloring is the same as that in Figure|4] 



pressure for = 8/3 is the only solution that goes to zero at 
a finite value, thereby forming a boundary layer of thickness 
^ ~ 1, or ~ 1/r. Above r; = 8/3 the pressure profiles 
begin to decrease more gradually again, and for 77 ^ 4 the 
profiles again approach a constant. 

Physically, this has very interesting implications. For rj = 
8/3, the pressure dropping to zero at a fixed value of ^ implies 
that all of the jet material has piled up in the outer region of 
the boundary layer. Thus the jet is a hollow cone: the ma- 
terial is concentrated against the outer wall of the jet, and 
the inner region is essentially evacuated. This "edge pileup" 
has been studied as a possibility in quasar jets previously, e.g. 
IZakamska et alj r2008. Observationally, Zakamska et al. sug- 
gested this phenomenon as an alternative explanation for ob- 
served edge-brightening in jets, which is more commonly in- 
terpreted as a result of the Kelvin-Helmholtz instability occur- 
ring at the jet boundary. 

For values of 77 near 8/3, while the material is not im- 
mediately all piled up against the wall, the structure is not 
dissimilar. The region against the outer wall has the highest 
pressure, and the sharp pressure decrease interior to the outer 
wall suggests that most of the material is still in the outer re- 
gion. While not truly a hollow cone, this structure could still 
account for observed edge-brightening in jets. Moving in ei- 
ther direction away from 77 = 8/3, however, this region be- 
comes broader and the jet becomes more evenly distributed. 

For large and small enough values of rj, the pressure gra- 
dient becomes small and the deviation of the pressure within 
the shock from the external pressure is minimal, suggesting 
that the Kompaneets approximation is a valid approximation 
in this region. For values of 77 nearer to 8/3, however, the 
Kompaneets approximation is clearly not a reasonable one. 

The fact that the pressure profile begins level at 77 = 2, 
steepens as 77 approaches 8/3, and then levels out again as 77 
approaches 4 can be explained as the result of a tradeoff be- 
tween two opposing effects: as 77 increases, the degree of col- 
limation decreases. This causes a decrease in the centripetal 
force, and therefore the pressure gradient should decrease ac- 
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cordingly. But an increase in r; also means an increase in the 
rate of jet acceleration (as measured by the increase of F with 
r), which will result in a steeper pressure gradient across the 
boundary. This acceleration effect appears to be stronger in 
the region where 2 < 77 < 8/3, but the coUimation effect wins 
out in the region where 8/3 < < 4. The two effects are in 
balance when ri — 8/3, leading to maximal compression of 
the gas against the wall of the jet. 

This dividing behavior at r; = 8/3 is not entirely unex- 
pected, as can be shown by a quick calculation. The power 
carried by a pressure-dominated jet is given by 

L oc pF'^A, (18) 

where p is the pressure and A is the cross-sectional area of 
the region carrying most of the power If we assume that 
the power is concentrated in an outer layer of width — 
1 /r, then viewed end-on the cross-sectional area of that ring 
would be A oc r^/F, supposing that the jet is conical to first 
order with very slow collimation. Using the external pressure 
profile and employing the relativistic Bernoulli equation, we 
see then that 

L oc r-^/^A = r-^''/*+2. (19) 

If we demand that the jet power be distance-independent, 
thus containing a finite amount of energy in the boundary 
layer, then this scaling implies that 77 = 8/3 must be true. 

For 77 < 8/3, we have L scaling with radius to some pos- 
itive power, such that for r — >• 00, L — >• 00. Thus, given a 
roughly conical jet, the power carried cannot remain constant 
with 77 < 8/3. Reexamining Eq il9h however, we can see that 
if the area scaled as some smaller power of r rather than as 
described, then the jet power could be maintained as a con- 
stant despite the lower value of 77. Thus, for 77 < 8/3, the 
jet must become more strongly collimated, causing the cross- 
sectional area to grow more slowly than in the conical case, 
in order for the jet power not to diverge. 

For 77 > 8/3, it would appear that i — > as r — > cx). 
This is misleading, however, given that for this range in 77, the 
integrals over ^ of the energy contained within the boundary 
layer diverge. Thus, in a total spatial integral of L over both ^ 
and r, the r-scaling of L — > and the ^-scaling of L — >• 00 can 
combine to counteract each other, making it possible to obtain 
physical solutions under which a finite amount of energy is 
contained within the boundary layer after all. 

3.2 Extension to the Kompaneets Model 

This self-similar construction provides an important view of 
the behavior within the boundary layer, but to properly un- 
derstand observations of the large-scale AGN jets that we seek 
to model, we must examine how the boundary-layer physics 
fits into that of the jet as a whole. 

To this end, we now attempt to extend the Kompaneets 
model by repeating §2 with the addition of a pressure gradi- 
ent across the boundary layer The form of this added pres- 
sure gradient will come from our self-similar boundary-layer 
model, and in this section we work in cylindrical coordinates 
to facilitate the matching. 

Applying the self-similar model to the structure devel- 
oped in §2 inherently renders the problem no longer self- 
similar, as length scales are introduced into the problem. 
Nonetheless, the combination of methods is instructive in the 
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Figure 7. The shock front and contact discontinuity for two cases: one 
in which the jet originally closes (a) and one in which it originally re- 
mains open (b) . Left plots demonstrate the jet under the Kompaneets 
approximation, from §2. Right plots indicate the new position of the 
shock front after the addition of the pressure gradient. The contact 
discontinuity remains unchanged in this approximation. 



description of general trends that are expected when adding 
a pressure gradient to the Kompaneets model. 

Because ^ is a function of distance from the contact dis- 
continuity {9c — 6), using the self-similar pressure solutions in 
the model in §2 requires that we already know the shape of 
this outer wall. As the self-similar model was developed with 
the assumption of an approximately conical contact discon- 
tinuity, we solve this problem by fitting a conical solution to 
each contact discontinuity found in §2 using the Kompaneets 
approximation and then using this fit as a fixed input for the 
location of the outer wall. As was shown in §2, approximating 
the contact discontinuity as conical is justified in the case of 
jets where the shock never closes to the axis, or in the region 
of the jet before the shock closes. It is this regime that we 
study. 

We examine three representative cases — one for r) < 
8/3, one for 7; = 8/3, and one for 77 > 8/3 — and use these to 
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infer the general behavior of the jet shape after the addition 
of a pressure gradient. 

For r; — 7/3, the physical self-similar solution is that for 
which /i = 4.575. From this we fit an anal)d:ic function to the 
numerical solution for g, of the form 

'^^^ ^ "+(?T7) + (cT^' '''' 

where a, b, c, d, and e are constant parameters determined by 
the fit. We then modify the expression for the pressure within 
the boundary layer so that ps oc r^'' g'^{^). We find that the 
results display the expected physical effect of adding this pres- 
sure gradient: the shock front moves outward to match the 
decreased pressure. 

For = 8/3, the solution for was analjrtic: 

9(0 = (21) 

Using this function to add a pressure gradient into the 
case from §2, we again find that the shock wall moves out- 
ward to match the decreased pressure function. Two exam- 
ples are shown in Figure [7| one in which the shock front 
initially closed to the axis in the Kompaneets approximation, 
and one in which it initially remained open. Due to the ap- 
proximation of a conical contact discontinuity breaking down 
beyond the point where the shock front initially closes, we 
cannot draw conclusions from these models in that regime. 
Nonetheless, it is clear that the addition of a pressure gradi- 
ent either causes the shock front to converge to the axis at 
a distance further down the jet axis, or else it prevents the 
shock front from ever reaching the axis. Either way, one can 
see from Figure [7] that the addition of the pressure gradient 
causes the shocked boundary layer to become thinner. 

For = 3, we can again fit the function given by Eq I I20D 
to the numerical solution of g(^). This time the physical so- 
lution corresponds to /i = 1/2. Again modifying the Kompa- 
neets solution of §2 with this pressure gradient, we find that 
here too the shock moves outward and forms a thin bound- 
ary layer. As with the 77 = 8/3 solutions, cases that originally 
closed to the axis in the Kompaneets approximation either 
close further down the axis or no longer close. 

We would hope that the energy within the shocked layer 
increases with radius, despite the message of the luminosity 
scaling worked out at the end of §3. As a test of this, for ri — 3, 
we can integrate the total energy contained in the boundary 
layer in a slice across the jet at a fixed z: 

Lex / pT^RdR. (22) 

Comparing this enclosed energy at a few different values 
of z (see Figure|8]l, we can see that although the total energy 
within the layer would be decreasing if ^ were fixed along the 
shock (the dashed curves show the fixed-^ energy through 
a given value of z), the net energy instead increases as we 
go to higher z because we move to different values of ^ in 
the process of traveling along the shock. Thus the solution is 
entirely physical: the energy behind the shock increases with 
an increase in distance from the source, as is required. 




Figure 8. The total integrated energy contained within the boundary 
layer at a given height 2 for r; = 3. L is given in arbitrary units. The 
dashed curves through the data points correspond to what the total- 
energy curve through each data point would be if g were instead a 
fixed value along the shock. 

4 CONCLUSION 

We have evaluated the shape of the shocked boundary layer of 
a hot jet with an ultrarelativistic equation of state, in the case 
that it is injected into an ambient medium that has a power- 
law pressure profile of p oc r^^ with 2 < 77 < 4. Using the 
shock jump conditions and momentum conservation and as- 
suming that the pressure remains constant across the bound- 
ary layer, we found that whether the shocked layer closes to 
the axis — and where it closes, if it does — is dependent upon 
the values of the initial pressure ratio, the power-law index, 
the initial bulk Lorentz factor, and the product of the initial 
opening angle and initial Lorentz factor. 

We also found that, in the Kompaneets approximation, 
the contact discontinuity asymptotes to the shape Rc oc 2;'"'* 
in the case where the shock has closed to the axis, but to a 
conical shape, Rc oc z, in the case of a non-closing shock or 
in the region of the jet before the shock has closed. 

Due to the expectation of a pressure gradient arising from 
the centripetal force created by the slight collimation of the 
jet, we then created a self-similar model of the boundary layer 
that allows for a pressure gradient across the layer. From this 
we found solutions for which pressure decreases monotoni- 
cally across the boundary layer for all values of the power-law 
index r/. 

An 77 of 8/3 constitutes a special solution where the pres- 
sure goes to zero at a finite distance from the outer wall. For 
values of rj greater or less than 8/3, the pressure approaches 
zero asymptotically. The curve for pressure steepens for val- 
ues of rj near 8/3 and approaches a constant value for ry — >• 2 
or ?7 — >■ 4. 

The drastic decrease of the pressure inwards from the 
contact discontinuity for values of 77 near 8/3 suggests that 
most of the material is pushed up against the outer wall, cre- 
ating a hollow-cone structure for the jet that may be observ- 
able as edge-brightening. It also indicates that the Kompa- 
neets approximation is not valid for these values of rj. For 
values of 77 near 2 or 4, however, the Kompaneets approxima- 
tion may be reasonable, since the pressure decrease is very 
gradual. 

To better understand these results, we then revised the 
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Kompaneets solutions to include the pressure gradient of the 
self-similar boundary-layer model. We held the outer wall 
fixed and solved for the new position of the shock front by 
pressure-matching across the shock jump. We found that the 
addition of the pressure gradient caused the shock front to 
move outward in order to match the lower pressure, resulting 
in a thinner boundary layer and preventing the shock from 
closing to the axis in places where it originally had done so. 
We also confirmed that the total energy contained within the 
boundary layer was an increasing function of height z, as is 
required. 

The inherent difficulty of examining this problem ana- 
lytically required us to make several simplifying assumptions 
over the course of this work. One such assumption, made 
in §3 to obtain the boundary-layer solutions, is that all the 
material within the boundary layer is on the same adiabat. 
This would be true if the majority of the material entered the 
boundary layer at roughly the same location, however this 
is not the case in general. But for the boundary-layer solu- 
tions where the pressure drops rapidly (77 near 8/3), most 
of the material must be pushed up against the contact dis- 
continuity early on. This suggests that it must all have en- 
tered the boundary layer near the jet base, which demon- 
strates self-consistency with the same-adiabat assumption for 
this regime. Nonetheless, in future work we intend to explore 
the possibility of treating the jet material on different adia- 
bats. 

It is widely believed that magnetic effect s contribute 
to th e coUimation of relativistic jets (see e.g. iKomissarovl 
Il999h . While the models discussed here are purely hydro- 
dynamic, they can have important applications both in in- 
terpreting the data obtained in three-dimensional general 
relativistic magnetohydrodynamics simulations that provide 
a self-co nsistent description of the je t launching mecha- 
nism Ce.g. fBeckwith et"aLll2008 . 2009; M cKinnev & BlandfordI 
12009") and in improving the physical content of boundary con- 
ditions imposed in simulations that study th e role of magnetic 
fields in large-scale jet coUima tion (see e.g. iKomissarov et alj 
I2OO7I l2009t lKomissarovll201 ll) . 

In the former case, the Poynti ng-flux dominated jet i s 
sheathed by an unbound outflow jHawlev & KrolikI [20061) . 
which may have important implications for the operation 
of current-driven instabilities ( McKinnev & Blandford 2009). 
However, the numerical resolution used in these simulations 
is generally insufficient to adequately resolve the sharp den- 
sity and pressure gradients present in these regions. Combin- 
ing the models presented here with the conditions present 
in the ambient medium at the base of the jet (e.g. the disk 
corona) in these simulations will allow us to assess the extent 
to which numerical resolution affects simulation outcomes. 

In the case of simulations that study the role of magnetic 
fields in large-scale jet coUimation, the work presented here 
will allow a more complete physical treatme nt of the effects of 
the ambient medium. In these simulations dKomissarov et alj 
I2OO7', '2OO9I; IKomissarovl I2OI I), this boundary is generally 
treated as a rigid wall, where fluid quantities (such as gas 
density and pressure) are simply copied across the boundary, 
while the normal component of vector quantities (e.g. velocity 
and magnetic field) is reflected. The models presented here 
will allow improvement of this treatment by specifying jump 
conditions for hydrodynamic quantities at this boundary, con- 



sistent with the shape of the wall — thereby including the 
collimating effect of the ambient medium. 

These calculations mark the first step towards a more 
complete treatment of both the action of the external medium 
and magnetic effects in collimating and accelerating relativis- 
tic jets. The next step in this work is to include the effect of 
a toroidal magnetic field, allowing for magnetic confinement. 
Together, this work and its magnetized extension will provide 
us with an equilibrium model from which we can explore in- 
stabilities and radiative mechanisms of the jet near its base, 
leading to improved constraints on jet dissipation and associ- 
ated radiative signatures. Through this work we hope to fur- 
ther our understanding of jet coUimation and, in particular, 
boundary-layer behavior. 
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